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Hybrids of graphene and two dimensional transition metal dichalcogenides (TMDC) have the 
potential to bring graphene spintronics to the next level. As we show here by performing hrst- 
principles calculations of graphene on monolayer M0S2, there are several advantages of such hybrids 
over pristine graphene. First, Dirac electrons in graphene exhibit a giant global proximity spin-orbit 
coupling, without compromising the semimetallic character of the whole system at zero field. Re¬ 
markably, these spin-orbit effects can be very accurately described by a simple effective Hamiltonian. 
Second, the Fermi level can be tuned by a transverse electric field to cross the M0S2 conduction 
band, creating a system of coupled massive and massles electron gases. Both charge and spin trans¬ 
port in such systems should be unique. Finally, we propose to use graphene/TMDC structures as a 
platform for optospintronics, in particular for optical spin injection into graphene and for studying 
spin transfer between TMDC and graphene. 


Graphene spintronics [T] has relied exclusively on elec¬ 
trical spin injection [IS]. Combining graphene with 
semiconducting two-dimensional TDMC [5] can open 
new venues for spintronics applications 13 |7|. Indeed, 
TMDC are becoming increasingly popular in optoelec¬ 
tronics as sensitive photodetectors [8] or, forming lat¬ 
eral heterostructures [H |T0], as two-dimensional solar 
cells m- Important, TMDC have a sizeable spin-orbit 
coupling and lack space inversion symmetry. As a re¬ 
sult, their band structure [12] allows for a valley resolved 
optical spin excitation by circularly polarized light [I3f 
im. TMDC can thus facilitate optical spin injection into 
graphene, in hybrid structures. 

Efficient growth of M 0 S 2 on graphene has already been 
demonstrated [TBUTS] . It was reported that graphene 
on M 0 S 2 is ultrafiat, having large mean free paths [19]; 
angle-resolved photoemission found an intact Dirac point 
but a strong hybridization elsewhere in the tt system [20] . 
Technological potentials for these hybrid structures are 
already being discussed [21], mainly as a basis for non¬ 
volatile memory [22], sensitive photodetection [23], and 
gate-tunable persistent photoconductivity m- Recently, 
the spin Hall effect in graphene on few-layer WS 2 was ob¬ 
served at room temperature [25] . 

In this paper we establish by first-principles calcu¬ 
lations fundamental electronic properties and the spin- 
orbit fine structure of the graphene Dirac bands for 
graphene on monolayer M 0 S 2 , and introduce an effec¬ 
tive spin-orbit Hamiltonian which explains the proxim¬ 
ity induced spin splittings of the Dirac states. We show 
that the induced spin-orbit coupling is giant, being 20 
times more than in pristine graphene. We also discuss 
the field effect on the band offsets of the two materials. 
Finally, we present possible experimental schemes to per¬ 
form optical spin injection into graphene and study spin 
tunneling from TMDC through graphene. 


First-principles results. To establish the electronic 
and spin properties of graphene on M 0 S 2 we used first- 
principles methods based on density-functional theory 
[26] . see Methods. To reduce structural strain we con¬ 
structed a large supercell of 59 atoms, comprising a 3 x 3 
supercell of M 0 S 2 and a 4 x 4 supercell of graphene, with 
the residual lattice mismatch of 1.4%. The relaxed inter¬ 
layer distance between graphene and M 0 S 2 is 3.37 A, see 
Fig.[^). In this supercell the K point of M 0 S 2 is mapped 
to the F point in the reduced Brillouin zone. The calcu¬ 
lated electronic band structure is shown in Fig. [^). The 
Dirac cones of graphene are nicely preserved, with the 
projected Dirac point (which is also the Fermi level) be¬ 
ing slightly below the conduction band edge of M 0 S 2 . 
The closeness of the Dirac point to the conduction band 
of M 0 S 2 enhances screening, which can substantially in¬ 
crease the mean free path in the graphene layer, as re¬ 
cently shown experimentally m- 

The band offsets between graphene and M 0 S 2 can be 
controlled by an external electric field applied transverse 
to the layers. This is demonstrated by our first-principles 
calculations in Fig. |^), where we present Ac, the differ¬ 
ence between the conduction band minima of M 0 S 2 and 
graphene. At negative fields (pointing towards M 0 S 2 ) 
the offset increases, leaving both layers neutral. How¬ 
ever, positive fields shift the Dirac point above the con¬ 
duction band minimum of M 0 S 2 and populate graphene 
with holes and M 0 S 2 with electrons. The Fermi level 
crosses both the valence band of graphene and the con¬ 
duction band of M 0 S 2 , see Fig. IT) and c). This field 
effect can establish a unique system in which massless 
Dirac electrons are coupled with a conventional 2d elec¬ 
tron gas m- 

We now zoom in on the Dirac point at K to see how 
the electronic spectrum of graphene deforms in the pres¬ 
ence of M 0 S 2 . This fine structure is shown in the inset to 
Fig. [^). There are two important effects: First, an or- 
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FIG. 1 : Calculated electronic and spin properties of graphene 
on two-dimensional M0S2. (a) The supercell used in first- 

principles calculations, (b) Calculated band structure along 
high symmetry lines. The inset is a zoom to the hne structure 
of the low energy bands at the Fermi level, around the Dirac 
point. Bands with positive (negative) z component of the spin 
are shown in red (blue). The sublattice character {A and B) 
is also indicated, (c) Spin textures for the four bands of the 
inset in b). 


bital band gap opens, due to the breaking of the graphene 
pseudospin symmetry. On average, atoms A and B in the 
graphene supercell see a different environment coming 
from the M 0 S 2 layer. This orbital gap is there even in the 
absence of spin-orbit coupling. It arises from the effective 
staggered potential induced by the pseudospin symmetry 
breaking. Second, spin-orbit coupling combined with the 
broken space inversion symmetry lifts the spin degener¬ 
acy of the Dirac valence and conduction bands and leads 
to the appearance of four distinct bands. This splitting 
is on the meV scale, which is giant when compared to 
the 24 /ieV spin-orbit splitting in pristine graphene [28]. 
The inset also shows the orbital character of the bands at 
K: while the valence states are formed at the B sublat¬ 
tice, the conduction states live on A. The same orbital 
ordering is at K'. 

Another important characteric of the Dirac states is 
their spin texture. This is plotted in Fig. El) for the 
four bands from the inset of Fig. [^). Directly at K the 
spins are pointing out of the graphene plane, alternating 
up and down. Increasing the momentum away from iF, 
the spins acquire a winding in-plane component, either 
clockwise or counterclockwise, suggestive of the strong 
Rashba effect. At K' the spins are reversed. 

Effective Hamiltonian. Can we understand these 
proximity-induced changes in graphene’s band structure 



FIG. 2 : Field effect on graphene/MoS2. (a), Calculated offset 
Ac from the conduction band minimum of M0S2 to graphene, 
as a function of an applied transverse electric held. At positive 
helds electrons are transferred from graphene to M0S2, estab¬ 
lishing a massless-massive electrons bilayer, (b) Net Fermi 
energy for graphene, as the difference of the valence band 
maximum Egv of graphene and the system Fermi level, Ef- 
At positive helds the valence band of graphene becomes pop¬ 
ulated. (c) Net Fermi energy for M0S2, as the difference be¬ 
tween the system Fermi level, Ef, and the conduction band 
minimum, Ec of M0S2. The conduction band of M0S2 be¬ 
comes populated at positive held, rehecting the population of 
holes in graphene in (b) as the whole system is neutral. 


from an effective model? The answer is not obvious since 
not only sublattices A and B differ, but even the sites 
that belong to the same sublattice see different local en¬ 
vironments in the supercell. Surprisingly, an effective 
symmetry-based Hamiltonian with graphene orbitals in 
the presence of pseudospin inversion symmetry breaking 
gives a remarkably good description. The model builds 
on the orbital Hamiltonian for pristine graphene which, 
close to K{K') points, is 


Ho = hv^ii^cr^kx + CTyky). ( 1 ) 

Here vf is the Fermi velocity of graphene, kx and ky 
are the Cartesian components of the electron wave vector 
measured from iF(iF'), parameter n = 1 (—1) for K {K')^ 
and ax and ay are the pseudospin Pauli matrices acting 
on the two-dimensional vector space formed by the two 
triangular sublattices of graphene. Hamiltonian Ho de- 
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scribes gapless Dirac states with the conical dispersion 
eo = uhvY\k\ near the Dirac points; u = 1(—1) for the 
conduction (valence) band. 

The staggered potential describing the effective orbital 
energy difference on A and B sublattices of graphene on 
M 0 S 2 enters via the Hamiltonian, 

■Ha = A(Tj;So, (2) 

where is the pseudospin Pauli matrix and sq is the 
unit spin matrix; A is the proximity induced gap of the 
Dirac spectrum. Another consequence of the pseudospin 
inversion asymmetry is the sublattice-resolved intrinsic 
spin-orbit coupling. Indeed, the intrinsic coupling acts 
solely on a given sublattice: it is a next-nearest neighbor 
hopping [29]. We describe it in our model with param¬ 
eters and Xf for sublattices A and H, respectively. 
The corresponding proximity induced spin-orbit coupling 
Hamiltonian close to K{K')^ 

'Hso = Af-[((T2 + cro)/2]KS2 + Af [(cr^ - ao)/2]Ks^, (3) 

is a generalization of the McClure-Yafet Hamiltonian for 
graphene mm- We denote by Sz the spin Pauli matrix, 
while by ao the unit matrix acting on the pseudospin 
(sublattice) space. If = Af, the main effect of the 
intrinsic spin-orbit coupling is to enhance the anticross¬ 
ing of the pristine graphene Dirac cones [28] , leaving the 


spin degeneracy intact. However, if A^ 7 ^ Af, as in our 
case, the spin degeneracy gets lifted by this intrinsic term 
already, reflecting the loss of space inversion symmetry. 

Placing graphene on M 0 S 2 also breaks the lateral mir¬ 
ror symmetry, giving rise to the Rashba type spin-orbit 
coupling [29] . 


■Hr = XniKa^Sy - cJyS^), (4) 


where Ar is the Rashba parameter and Sx^ Sy are the 
spin Pauli matrices. In the hopping language, the Rashba 
coupling is the nearest-neighbor spin-flip hopping, con¬ 
tributing further to the spin splitting of the bands, and 
defining the spin quantization axis for each Bloch state, 
away from the time reversal points P and M. 

Hamiltonian Ho + Ha + ^so + fully describes 
graphene’s bands at K{K'). Its eigenenergies are 




1 + Z^/i 


uA 




1 — Ujii 


Af + Af - vJ{2A - Af + Af )2 + 16A| 


, ( 5 ) 


where /i = 1(—1) for spin up (down) branches. The ex¬ 
pectation values of the spin along z for the corresponding 
states are given by 


IJ.hvh 

\Sz)vil = r) 


1 ^ Ufl 1 — Ufi 


Using the formulas for the eigenenergies, Eq.([^, and 
for the spin expectation values, Eq.(|^, we can alge¬ 
braically extract the orbital band gap A and the three 
spin-orbit parameters A^, Af, and Ar by comparing to 
our first-principles data for the fine structure at H, see 
the inset to Eig. Et)- The extracted parameters are 
shown in Eig. [^ as a function of the applied transverse 
electric field. The orbital proximity gap A is about 0.5 
meV in zero field. In fields greater than 0.5 V/nm, 
the gap exhibits a steep increase, which is related to 
the transfer of the electronic charge from graphene to 
M 0 S 2 , see Eig.[^ The proximity spin-orbit parameters in 
Eig.[^) and c) are about 0.2 meV, which is 20 times more 
than in pristine graphene [28]. Similar giant values of 
spin-orbit coupling in graphene are induced by hydrogen 
adatoms [3TII33] . and even more by fluorine [34], though 
the mechanisms are different. Unlike in the adatom cases 
in which the induced spin-orbit coupling is only local, in 
our case the giant coupling is global. While the intrinsic 


2A - Af - Af 
\J (2A — Af + '^p)^ T IbA^ 


(6) 


parameters Af and Af change rather moderately with 
applying the electric field, the Rashba parameter Ar, see 
Eig. 1 ^), more than doubles in increasing the field from 
-2 to 2 V/nm. 

What is the origin of the induced giant spin-orbit cou¬ 
pling in graphene on M 0 S 2 ? We trace the enhancement 
to the hybridization of the carbon orbitals with the d- 
orbitals of Mo. We find only 0.3% of d-orbitals at the K 
point by analyzing the calculated density of states. But 
when we turn the spin-orbit coupling on Mo atoms in the 
supercell off, the orbital gap in zero field remains almost 
unchanged (A = 0.506 meV), while the spin-orbit param¬ 
eters drop to their pristine graphene values Af = 24 /ieV, 
Af = 23 /ieV and Ar = 10 /ieV, which are, curiosly, also 
determined by d orbitals, but from carbon atoms [28] . 

Away from H(H'), the spin splittings depend on the 
momentum. In order to describe our first-principles data, 
we add the PIA (pseudospin inversion asymmetry) spin- 
orbit coupling term [32] which, like the intrinsic coupling. 
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FIG. 3 : Calculated effective Hamiltonian parameters at the K 
point: (a) hybridization gap A, (b) sublattice resolved intrin¬ 
sic spin-orbit coupling and Af (they have opposite signs), 
and (c) Rashba parameter Ar, as functions of the applied 
transverse electric field. 


represents the next nearest neighbor hopping, but with a 
spin flip. The full model Hamiltonian describes the data 
perfectly, see Supplementary information. 

Although our effective model should capture the basic 
physics of graphene on TMDC, the extracted parame¬ 
ters are for the specific supercell of graphene on M 0 S 2 . 
Certainly, taking an even larger cell that could further 
reduce strain, or twisting the two layers as could hap¬ 
pen in experiments, would lead to a different set of pa¬ 
rameters, although the orders of magnitudes would likely 
stay. In a macroscopic experimental structure we ex¬ 
pect Moire patterns which would transform our Hamil¬ 
tonian into a Hamiltonian density, with an orbital gap 
and spin-orbit fields, perhaps even averaging some of the 
parameters (such as A and the difference between Af 
and Af) to zero. Our extraxted parameters can then 
be viewed as effective standard deviations of the spatial 
variations, suitable as input for charge and spin transport 
model calculations for such samples. We also expect that 
graphene on TMDC could produce superlattice features 
as in graphene on hBN [35] . 

Optospintronics. We propose graphene-TMDC hy¬ 
brids, such as the one studied above based on M 0 S 2 , as an 
ideal platform for optospintronics. In Fig.[^) we give an 
optical spin injection scheme into graphene. A circularly 
polarized light, tuned to the band gap of TMDC, excites 
electron spins by optical orientation [SlISS]- In effect, the 
light produces spin-polarized excitons which dissociate 
into spin-polarized electrons and holes. As in the recent 
optical experiment [24], we expect that electrons will be 
transfered to graphene, leaving holes behind in TMDC, 
although in which way electrons and holes split may de¬ 
pend on the TMDC material as well as on gating. The 



FIG. 4 : Optospintronic schemes for graphene-TMDG hybrids, 
(a) Optical spin injection into graphene, facilitated by the 
semiconducting TMDG. A circularly polarized light excites 
spin-polarized electrons in the semiconductor. The spin is 
transferred to graphene where it can be detected as a Hanle 
signal by the ferromagnetic electrode, (b) Spin transfer be¬ 
tween two different TMDGs, encapsulating graphene. Gircu- 
larly polarized light tuned to the band gap of the top material 
excites electron spins which can tunnel to the lower material, 
exhibiting a circular luminescence peaked at its band gap fre¬ 
quency. 


spin-polarized electrons (or holes) diffuse in graphene. 
One can detect this spin accumulation either optically, 
by observing a circular polarization of the photolumi¬ 
nescence [36] elsewhere in graphene on TMDC, or elec¬ 
trically. The latter is illustrated in Fig. |^): a ferromag¬ 
netic electrode on top of graphene detects the presence of 
the spin accumulation in graphene mm- Spin precession 
in graphene can be observed as the Hanle signal (which 
is not possible to see in the spin-valley coupled TMDC 
EU), by applying an external magnetic field transverse 
to the injected spin, providing Larmor precession [7j. 

Spin transport per se in graphene-TMDC bilayers 
should be fascinating. The presence of the giant, effectiv- 
elly uniform spin-orbit fields should give large spin Hall 
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signals, even greater than in hydrogenated graphene [33] . 
Most important, as our calculations show, the spin, like 
charge, properties of these structures are expected to be 
highly field tunable. The fascinating prospect of realiz¬ 
ing the massive-massless electron gas coupling of the two 
electron gases, if the Fermi level is positioned in both 
band structures, calls for new theories of spin transport 
and spin relaxation in such hybrid systems. 

To demonstrate spin tunneling from a TMDC through 
graphene one could use a sandwich structure, as pictured 
in Fig. [^). The two semiconductors have different 
band gaps, allowing to discern the photoluminesence sig¬ 
nals from the top and bottom layers. If the spin pump¬ 
ing light is tuned to the band gap of the top layer, the 
spin-polarized carriers would be excited there and tun¬ 
nel through graphene to the bottom layer, in which they 
would recombine and emit circularly-polarized light with 
the frequency characteristics of the bottom material. One 
can envision influencing the signal with a transverse mag¬ 
netic field, allowing for a Hanle effect Another possibility 
is to measure the accumulated spin in the bottom layer 
using the magneto-optical Kerr effect, as in recent exper¬ 
iments on monolayer M 0 S 2 [38] . 

Conclusion. We have established by first-principles 
calculations a strong effect of M 0 S 2 on the spin properties 
of graphene, predicting a giant and field-tunable proxim¬ 
ity spin-orbit coupling for Dirac electrons. We have in¬ 
troduced an effective spin-orbit Hamiltonian to describe 
the electronic states around the Fermi level, fitting per¬ 
fectly the first-principles data. We have also showed that 
gating can tune the band offsets of the two layers, al¬ 
lowing to realize the unique system of coupled massless 
and massive electron gases. Finally, we have proposed 
to use graphene on TMDC as a platform for optospin- 
tronics with graphene-based two-dimensional materials 
structures. 

We thank T. Korn, C. Schuller, C. Stampfer, B. 
Beschoten, T. Muller, and B. Ozyilmaz for use¬ 
ful descussions and hints regarding possible experi¬ 
mental realizations of optospintronics with graphene- 
TMDC structures. This work was supported by 
DFG SFB 689, GRK 1570, and by the EU Sev¬ 
enth Framework Programme under Grant Agree¬ 
ment No. 604391 Graphene Flagship. 

Supplementary information 
Methods 

Structural relaxation and electronic structure calcula¬ 
tions were performed with Quantum ESPRESSO [39] . 
using norm conserving pseudopotentials with kinetic en¬ 
ergy cutoff of 60 Ry for wavefunctions. Eor the exchange- 
correlation potential we used generalized gradient ap¬ 
proximation m- The supercell containing a 3 x 3 super¬ 


cell of M 0 S 2 and 4x4 supercell of graphene was embeded 
in a slab geometry with vacuum of about 13 A, with a 
dipole correction im, which is crucial to get accurate 
band offsets between the Dirac point and the conduc¬ 
tion band minimum of M 0 S 2 . The resulting structure 
has a lattice mismatch of 2.8% which we split equally be¬ 
tween graphene and M 0 S 2 by compressing graphene and 
stretching M 0 S 2 by 1.4%. The supercell has 59 atoms. 
The reduced Brillouin zone was sampled with 12 x 12 
k points. Atomic positions were relaxed using the quasi¬ 
newton algorithm based on the trust radius procedure in¬ 
cluding the van der Waals interaction which was treated 
within a semiempirical approach [42j 03] . The calculated 
work function on the graphene side is 4.12 eV, while on 
M 0 S 2 it is 4.41 eV. 

Spin eigenstates 

The normalized eigenstates at the K point of the 
Hamiltonian discussed in the manuscript in the basis a G 
s = IA t, A I, H t", ^ i) read 


t/)+_ = 

|0, 

4ARa+, 

0) 


V’++ = 

|1, 0, 

0, 

0) 


Ip — = 

o 

o 

0, 

1) 

= \Bi), 

^-+ = 

0, ia-Q-, 

4ArQ:_, 

0) 


where a± = 

1/+ ( 

3^ and 

-H 

= 2A + A; 

(2^ + 

Af + Af )2 + 

16A|. The 

eigenvectors 


ordered with increasing energy (using the extracted pa¬ 
rameters): _> 5 ++ >5 _> £_Analyzing the 

above eigenvectors we see that the valence bands are lo¬ 
calized on sublattice H, while the conduction bands on 
sublattice A. The z component of the spin alternates 
from band to band. This behavior matches the first- 
principles results, see inset to Eig.[^) in the manuscript. 

The top valence ip _and bottom conduction (/?++ states 

are pure pseudospin and spin states. On the other hand, 
spin-orbit coupling mixes spin and pseudospin of the out¬ 
ermost states. The eigenstates at K' have the same form, 
but opposite spins. 

PIA coupling: spin splitting away from K(K') 

To describe the calculated spin splittings away from 
A, we employ the PIA (short for pseudospin inversion 
asymmetry) spin-orbit term, introduced to study the ef¬ 
fects of spin-orbit coupling in graphene due to hydrogen 
adatoms [32] : 

^PIA ~ ^y^x) • (8) 

Here Apj^ and Apj^ are the spin-orbit parameters rep¬ 
resenting the average, Apj^, and differential, Apj^, PIA 
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coupling between the A and B sublattices. Like intrinsic 
spin-orbit coupling, PIA can be also represented by next- 
nearest-neighbor (same sublattice) hopping, but with a 
spin flip. The PIA terms turn the spin quantization axes 
of the electron states towards the graphene plane and 
add to the Rashba term for momenta away from the K 
point. 




electric field [V/nm] 


FIG. 5: Spin splitting away from K: emergence of PIA spin- 
orbit coupling in graphene on M 0 S 2 . (a) Spin splitting of 

conduction and valence bands from K (k = 0) in the direc¬ 
tion towards F, at zero electric field. Solid lines are model 
fits, symbols are first-principles results, (b) Fitted PIA pa¬ 
rameters at different transverse electric fields. 


To obtain Apj^ and Apj^ we add Hpia to the Hamil¬ 
tonian Bo + Ba + Bso + Hr, which is used directly at 
K in the manuscript, and fit to the first-principles data, 
keeping all other parameters as determined directly at 
K. In Fig. we plot the calculated spin splittings of 
the valence and conduction bands. The full model, with 
PIA, fits the first-principles data perfectly. The fits are 
Apj^ = —4.24 meVA, and Apj^ = 1.25 meVA. In Fig. 
we plot the two PIA parameters as functions of the ap¬ 
plied transverse electric field. Their tunability is enor¬ 
mous, much more than that of the Rashba hopping shown 
in the manuscript. As a final check we calculate, using 
the extracted parameters, the z components of the spin 
expectation values {s^) in the vicinity of the K point 
for the low energy states. The model is fully consistent 
with the first-principles data, as seen in the comparison 
plotted in Fig. 
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